Self- Assembled Triply Periodic Minimal Surfaces as moulds for 

Photonic Band Gap Materials 

L. Martm-Moreno*, F. J. Garcia- Vidal I' and A. M. Somoza * 
Departamento de Fisica de la Materia Condensada, ICMA-Consejo Superior de Investigaciones Cientificas, Universidad de 

Zaragoza, Zaragoza 50015, Spain. 
^ Departamento de Fisica Teorica de la Materia Condensada, Universidad Autonoma de Madrid, Cantoblanco, 28049, 

Madrid, Spain. 

* Departamento de Fisica. Universidad de Murcia. Aptdo. 4021, Murcia E-30080, Spain. 

We propose systems with structures defined by self-assembled triply periodic minimal surfaces 
(STPMS) as candidates for photonic bandgap materials. To support our proposal we have calculated 
the photonic bands for different STPMS and we have found that, at least, the double diamond 
and gyroid structures present full photonic bandgaps. Given the great variety of systems which 
crystalize in these structures, the diversity of possible materials that form them and the range of 
lattice constants they present, the construction of photonic bandgap materials with gaps in the 
visible range may be presently within reach. 
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A three-dimensional (3D) photonic bandgap material 
(PBGM) is a periodic dielectric system with an absolute 
frequency gap for electromagnetic (EM) waves [Q. The 
opening of photonic gaps is a delicate balance between 
both refraction index contrast and topology of the under- 
lying lattice. Since its proposal Q [|| and first construc- 
tion in the microwave range 0|, PBGMs have attracted 
much attention due to their many remarkable applica- 
tions especially for gaps in the visible range. However, 
material engineering of 3D periodic systems structured 
at the relevant length scale for optical gaps (hundreds 
of nm) is facing technological problems Another 
promising strategy is mesoscopic self-assembly: in this 
way partial bandgaps have been reported in 3D periodic 
structures synthesized from submicron colloidal particles 
[0 Within this route, an alternative idea has appar- 
ently been overlooked: the use of self-assembled triply 
periodic minimal surfaces (STPMS) as natural moulds 
to construct optical PBGMs. STPMS are present in dif- 
ferent systems like the periodic phases found in diblock- 
copolymers Band lipid/water systems Recently, 
the predicted ||l 1| periodic phases in microemulsions have 
also been observed p^ ]. 

All the binary and ternary systems mentioned above 
have the common property of being self-organized at a 
supramolecular level, with morphologies mainly dictated 
by the shape of internal interfaces separating different 
domains. In AB diblock-copolymers these internal in- 
terfaces are formed by the chemical bonds between the 
A and B blocks, whereas in lipid/water systems or in 
microemulsions (of for example water and oil) the am- 
phiphilic molecules form a monolayer interface that fully 
separates the polar compounds (water and amphiphilic 
heads) from covalent ones (oil and amphiphilic tails). 
From the physical point of view, the dominant factor that 
controls the shape of these internal interfaces is minimiza- 
tion of its surface free energy. This last condition is essen- 
tially equivalent to area-minimisation with fixed volume 
fractions and leads to the formation of constant mean 
curvature interfaces. When the different domains have 
approximately equal volume fractions, the mean curva- 
ture is zero. This basically explains why and when zero 



mean curvature interfaces are present in these systems. 
These kind of surfaces are known by mathematicians as 
minimal surfaces and, already in the last century [ p^ , 
it was found that some minimal surfaces can be joined 
together forming periodic structures in all three spatial 
directions. Subsequentently, it was discovered that 
minimal surfaces appear in a variety of real systems oc- 
curring in biology and materials science. 

Presently known systems with STPMS have three ap- 
pealing characteristics for an optical PBGM: i) Their lat- 
tice constants can be tuned by varying their constituents: 
real systems with STPMS have lattice constants ranging 
from one to hundreds of nm. ii) Their morphologies are 
bicontinuous with cubic symmetry, which is known [Q to 
be more appropriate to produce photonic bandgaps than 
the close-packed morphologies resulting in self-assembled 
colloidal crystals, iii) Flexibility in the composition, ei- 
ther by the appropriate choice of the constituents, or by 
modifications through selective chemical reactions, may 
provide the adequate refraction index contrast. 

However, although promising, these characteristics are 
not enough to secure the opening of photonic bandgaps: 
as previously mentioned, topology also plays a crucial 
role. So, in order to support our proposal of STPMS as 
possible PBGM, we have studied theoretically which ar- 
chitectures, if any, among the most commonly observed 
STPMS, develop bandgaps. We have performed calcula- 
tions for the structures defined by three different minimal 
surfaces: simple cubic (P), gyroid (G) and double dia- 
mond (D). Our computational scheme has two steps: (i) 
generation of the different structures and (ii) calculation 
of the photonic bands for each one of them. We have 
restricted ourselves to systems with two different com- 
ponents of equal volume fraction domains, characterized 
by their refraction indexes ni and n2 (for definiteness we 
take nil < 712), separated by a minimal surface. 

There are several mathematical techniques (Weier- 
strass parameterization, finite element methods, ...) p^ ] 
for the generation of triply periodic minimal surfaces. 
Based on physical grounds, it is expected that these pe- 
riodic surfaces could be generated via minimisation of a 
Landau-Ginzburg free energy functional llH], result that 
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has been proved rigorously Q . We have chosen this last 
approach, finding extremal configurations of the func- 
tional by solving numerically the Cahn-Hilliard equation 
This equation resembles the time evolution of the 
local concentration c(r) in the dynamical process of phase 
separation in binary mixtures of polymers and alloys. 
The local concentration is forced to verify a continuity 
equation, which guarantees fixed volume fractions during 
the full dynamical process. As the final result we obtain 
the stable (or metastable) time-evolved c{r) evaluated at 
a 3D regular grid of points covering the unit cell. 
Different triply periodic minimal surfaces, defined as the 
interface between both constituents, are obtained for dif- 
ferent initial configurations of c(r); known solutions are 
recovered by initializing c(r) near to the expected result. 
In Figs, la, lb and Ic we show, respectively, the P, G 
and D structures obtained via this method. 

As regards to the photonic band calculation, we have 
used a finite difference time domain (FDTD) method . 
This method requires a computational time which scales 
linearly with the number of mesh points inside the unit 
cell and, therefore, is specially suitable for systems with a 
dielectric constant that changes rapidly in space, like the 
ones we are studying. Basically, in the FDTD method an 
arbitrary initial EM wave with a wavenumber k is allowed 
to evolve both in time and space governed by a discrete 
version of Maxwell equations. After a certain run time, 
long enough to have good spectral resolution, the fre- 
quencies uj{k) of EM eigenmodes in the structure are ob- 
tained by Fourier transforming the time-propagated EM 
fields. For each structure we have performed calculations 
for values of N ranging from = 14 to = 40. We have 
checked that the results obtained converge rapidly with 
which gives us confidence on the quality and consis- 
tency of both structural and photonic calculations. 

It is worth pointing out that, although we are think- 
ing on applications in the optical regime, our results are 
valid for any lattice constant due to the scaling properties 
of the macroscopic Maxwell equations Gaps can be 
characterised by two quantities: the wavelength at the 
midgap frequency (Aq) in units of the lattice constant 
(a), times the smallest refraction index (ni), Ao/(?T.ia), 
and the "figure of merit" , /, defined as the quotient be- 
tween the frequency width of the gap and the midgap 
frequency. / is, at the same time, a measure of the gap, 
of the reduction of transmission of EM waves though a fi- 
nite number of layers, and of the robustness of the gap to 
small deviations from uniformity of the lattice constant 
throughout the crystal. 

We have calculated the photonic bands as a function 
of the index of refraction contrast n = n2/ni, for each 
one of the cited 3D periodic structures. For the P struc- 
ture we have not found photonic gap for any n. This 
result is in keeping with photonic band calculations for 
the, different but related, system of dielectric spheres as- 
sembled in a simple cubic lattice. Inset of Fig. la shows 
the low-frequency photonic bands for the P structure for 
a representative n — \/T3, showing the lack of gap in 
the spectrum due to some degeneracies in points of high 
symmetry. However, and this is the main outcome of our 
calculations, we do find gaps for the G and D structures, 
provided there is a minimum refraction index contrast be- 



tween the two materials. This can be clearly seen in the 
insets of Figs, lb and Ic where the low-frequency pho- 
tonic bands for the G and D structures are shown, for the 
same representative refraction index contrast n = vT3. 
Fig. 2a summarizes our results for the G structure, show- 
ing the dependence of both Ao/(nia) and / with n. Fig. 
2b shows the same quantities as Fig. 2a, this time for the 
D structure. Notice that, due to the fact that midgap 
wavelengths are approximately twice as large as lattice 
constants (an even larger for large n or large ni), pho- 
tonic bandgaps in the optical regime are expected for 
lattice constants of « 200 — 350 nm or even smaller. 

It is noteworthy that not only there is a minimum value 
of n for the appearance of photonic gaps, but also the 
quality factor / increases with n both for the G and D 
structures, thus making large refraction index contrast a 
desirable feature. The great flexibility in the manipula- 
tion of known systems with STPMS may provide routes 
to obtain the desired refraction index contrast. For ex- 
ample, in AB diblock copolymers, the optical proper- 
ties can be tailored modifying the A and B monomers, 
triblock copolymers providing even more flexibility. In 
lyotropic systems and microemulsions their liquid nature 
turns into an advantage. One of the liquids may be solid- 
ified (through freezing or polymerization) and the other 
one extracted. Moreover, the domains of one of the con- 
stituents may be used as a template for many different 
chemical reactions. For instance, in this way cubic struc- 
tures of silica with remarkable periodicity have been al- 
ready obtained The reported lattice constant was 
about 10 nm but it might be possible to increase this 
value, as pore sizes larger than 100 nm have been found 
in disordered silica mesostructures Extraction of 
the surfactant and substitution by another material with 
a large dielectric constant may lead to a PBGM in the 
ultraviolet regime. Even more interestingly, materials 
presenting STPMS have been used as hosts for synthe- 
sizing semiconductors with very large refraction index 
^l| . Based on all these, we believe that the construction 
of 3D PBGMs in the ultraviolet and optical ranges based 
on STPMS could be built using present day chemistry. 
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Figure Captions 

Figure 1. Three unit cells for the different triply peri- 
odic minimal surfaces analyzed in this paper: (a) simple 
cubic P, (b) gyroid G and (c) double diamond D. To the 
right of each structure, the corresponding photonic band 
structure is shown, for the representative case in which 
the two materials separated by the minimal surface have 
a refraction index contrast n = \/l3. The wave vector 
varies across the simple cubic Brillouin zone defined by 
the high-symmetry points: G(0,0,0), X(0,0,^), M(0,f ,f ) 
and R(^,^,^), a being the lattice constant of the system. 
EM frequency is in units of where — is the light ve- 

1 " nia ni " 

locity in the lowest refraction index material. Note the 
appearance of a complete photonic bandgap, highlighted 
in yellow, in the G and D band structures. 

Figure 2. Quality factor, / (full line) and midgap 
wavelength Aq (dashed line) as a function of the refrac- 
tion index contrast for the (a) gyroid and (b) double 
diamond structures. / is the quotient between the gap 
width and the midgap frequency. Aq is in units of the 
lattice constant a times the lowest refraction index n\. 
For both structures there is a minimum refraction index 
contrast to open a bandgap, this value being n ~ 2.5 for 
the G case and n « 2 for the D one. 
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This figure "figl.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/9810299vl 
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